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A set of fluid equations, taking into account the spin properties of the electrons and positrons in a magneto- 
plasma, are derived. The magnetohydrodynamic limit of the pair plasma is investigated. It is shown that the 
microscopic spin properties of the electrons and positrons can lead to interesting macroscopic and collective 
effects in strongly magnetized plasmas. In particular, it is found that new Alfvenic solitary structures, governed 
by a modified Korteweg-de Vries equation, are allowed in such plasmas. These solitary structures vanish if the 
quantum spin effects are neglected. Our results should be of relevance for astrophysical plasmas, e.g. in pulsar 
magnetospheres, as well as low-temperature laboratory plasmas. 
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I. INTRODUCTION 



II. GOVERNING EQUATIONS 



The magnetoplasma, first discussed by Alfven (T|], is today 
a mature research topic, with a wide range of applications. 
Recently, so called quantum plasmas in which the quantum 
properties of the plasma particles are taken into account, have 
received attention (see e.g. J2. S 0] for an up-to-date set of 
references). The collective motion of quantum particles in 
magnetic fields thus gives a natural extension to the classi- 
cal theory of magnetohydrodynamic s (MHD) in terms of so 
called quantum magnetoplasmas. The discussion of quantum 
plasmas has shown that quantum collective effects can have 
interesting consequences both in laboratory and astrophysical 
environments [! B S I! [M [III El El H E 
and part of the literature has indeed been influence d by exper- 
imental progress and techniques (see e.g. 11191 |20j |2lJ, |22j, l23l 
l24l l25ll . where also strong field effects are discussed). From 
the laboratory perspective, the motion of particles with quan- 
tum spin in strong fields, using e.g. intense lasers, has at- 
tracted interest as a probe of quantum physical phenomena 
llH Hzl |H |M [U 111. These studies are however mainly 
focused on single particle properties. In Refs. lHH how- 
ever, the kinetic properties of spin plasmas was investigated. 
Strong fields appear in pulsar and magnetar environments 
1I34I l35l I3r3l . Discussions of quantum plasmas in such envi- 
ronments can be found in Refs. II15i [m tutu. Moreover, 
studies taking both certain quantum electrodynamical effects, 
such as photon s plit t ing , as well as collective particle effects 
have been made IB7ll38ll . 



In the present paper, starting from the basic set of equations 
presented in Refs. [3] and [4], we derive the governing MHD 
equations for a pair plasma. These MHD equations takes into 
account quantum properties, such as spin, of the electrons and 
positrons. It is found that the nonlinear propagation of Alfven 
waves in such quantum MHD (QMHD) pair plasmas is gov- 
erned by a modified Korteweg-de Vries equation. This equa- 
tion is known to support solitary structures. Moreover, if the 
quantum spin effects are neglected, the solitary structures van- 
ish, making them true quantum solitons. The results should be 
of interest for both laboratory and astrophysical plasmas. 



In Refs. Jl and fl, it was shown that the evolution of spin 
plasma particles in a strongly magnetized environment is gov- 
erned by the equations for the particle densities tij and veloc- 
ities Vj 



and 
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where qj stands for the charge of the electrons (q e = — e) or 
positrons (q p = e), m = m e = m„ is the electron (positron) 
mass, ffj is the collisional contribution between species j and 
the second species, and Fgj is the total quantum force density 
(for a complete expression of this force density, see Refs. [3] 
and J3l). Here we will make use of the following expression 
for the quantum force on the electrons/positrons 
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where the first term is the gradient of the so called Bohm po- 
tential, the second term comes from the spin and B = |B|. We 
note that tanh(x) = Bin(x), where Bin is the Brillouin func- 
tion with argument 1/2 describing particles of spin 1/2. The 
temperature 7) is measured in units of energy. Furthermore, 
we have introduced the Bohr magneton jiig = eh /2m, where h 
is Planck's constant divided by 2n. 

The coupling between the quantum plasma species is me- 
diated by the electromagnetic field. The total magnetic field 
include both the classical contribution (from currents j = 
Y.qj n jV j) an d the s P m sources, such that Ampere's law reads 



VxB = ^i (j+jM) 
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(4) 



including the magnetization spin current jmj — V x 
(2qjrijjj.BS /h\qj\) for each species, where S is the spin vector. 
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Furthermore, we need Faraday's law 
V xE 



~dt' 



(5) 



We note that the spin current is determined by the spin vec- 
tor S. In general the spin vector follows a separate evolution 
equation, see Ref. (3|], but in the limit where the time-scales 
are much longer than the Larmor period, the spin vector can 
be approximated by 
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where B is a unit vector in the direction of the magnetic field. 
The magnetization spin current is thus given by 



njjiB tanh 
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III. 



ELECTRON-POSITRON PLASMA AND THE MHD 
LIMIT 



We now introduce the total mass density p = m(n, 
tip), the centre-of-mass fluid flow velocity V = m(n e V, 
ripVp)/p, and the current density j = — en e Y e 
assume that P e p = T e _ p n e p . Using these definitions, we im 
mediately obtain 



eripVp, and 



t + v.( P v) = o. 



(7) 



Assuming quasi-neutrality, i.e. n e « n p , we next add the 
momentum conservation equations for the electrons and 
positrons to obtain 
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where 11 = \(T e + T p )/2m e ] I + (m 2 /p)j ® j is the total pressure 
tensor in the centre-of-mass frame and 
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The collisional contributions have here cancelled due to mo- 
mentum conservation. Subtracting the momentum equations 
for the electrons and positrons, assuming ^ = t]en e j, where 
tj is the resistivity, we also obtain the generalized Ohm's law 
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where AT = T e — T p , and AFg = Fq e 
in the limit of quasi-neutrality, the non-spin quantum force 



F e p . We note that 



contributions to Ohm's law cancel, thus only leaving the spin 
quantum force, as long as the species have different tempera- 
tures. Eq. ( [Tot is valid during rather general conditions, and 
accounts for different temperatures of the species, as well as 
for short scale lengths of the order of the gyro-radius. From 
now on we will however limit ourselves to an equal temper- 
ature plasma, with scale lengths longer than the gyro-radius. 
Dividing by the density and taking the curl of Eq. ( [Tol l, we 
obtain 



= V x (V x B - rjj) 



(11) 



For the cases where the Alfven speed is much smaller than 
the speed of light, the displacement current in (01 can be ne- 
glected, and thus the current can be expressed in terms of the 
magnetic field. In this case, Eqs. (0, ([8]) and ( fTTT ) constitute 
a closed system for B, V and p . But in plasmas with strong 
magnetic fields, for example close to pulsars and magnetars, 
the Alfven speed approaches c, and the displacement current 
must be included in Ampere's law. However, we still obtain a 
closed system for the same equations (0, ([8]l and ( fTTT ), simply 
by eliminating the electric field using Ohm's law and writing 
the current as 



j = J_VxB-£o3-(VxB)-j M 
where from © we write the spin current as 
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for a quasineutral electron-positron plasma with equal tem- 
perature T . Note that in (fT2l > we have for simplicity omitted 
dissipative effects (77 — > 0). 



IV. ONE-DIMENSIONAL ALFVEN WAVES 

The spin terms have rather different properties than the 
standard MHD- terms. As a consequence, the spin force does 
not necessarily need to be as large as the ordinary j x B-force, 
in order to significantly affect the evolution of the system. 
In order to illustrate this specific property, we consider the 
weakly nonlinear evolution of shear Alfven waves propagat- 
ing at an angle to the magnetic field. In the absence of spin 
effects, the Alfven waves propagates with the Alfven speed 

1 /2 

ca = [c 2 Bq/ (c 2 /ioPo +5q)1 , where Bq is the magnitude of 
the unperturbed magnetic field, and po is the unperturbed den- 
sity. We consider the case with spatial dependence on a single 
coordinate % = xsin0 +zcos0. As a prerequisite, we first 
consider the linearized equations. Letting the unperturbed 
magnetic field lie along the z-axis, and combining Eqs. ([8]), 
(fTTT l and (12\ (with r\ — ► ) we, as a first approximation, 
obtain 
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Here c^.sp = c\/(l + 5 sp ) is the spin-modified 
Alfven velocity, with the spin modification deter- 
mined by 5 sp = h(of,tanh(iJ,BBo/T)/2mc 2 (Dc, where 

(O p = [(e 2 («op + n 0e) I 's^m] 1 ! 2 is the plasma frequency 
of the pair plasma, which comes from the part of the spin 
current proportional to (po/m)V x [tanh (ji B Bo/T) B]. Here 
we have introduced the cyclotron frequency co c = eBo/m 
ll45ll . For the sake of definiteness we also consider waves 
propagating in the positive ^-direction, such that we can 
use the expression d/dt = c& jSp cos Od/dE, in the linear 
nondispersive approximation. Furthermore, from now on 
we assume that the spin effects are small in the sense that 
<C 1, and also that the spins are weakly aligned with 
the external magnetic field such that we can approximate 
tanh{^ B Bo/T)^^ B B Q /T. 

As is wellknown, a weakly nonlinear evolution typically 
leads to wave steepening effects, and the necessity to include 
dispersion on an equal footing. Following Ref. 1I39I1 . but with 
the inclusion of spin terms, treating the Hall current, i.e. the 
first term of equation dTOI . as a small correction, weakly dis- 
persive effects are thus kept. Considering only the positive 
propagating waves, the linear wave operator for shear Alfven 
waves can then be generalized to 
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where small spin-contribution to the dispersive term in Eq. 
( Tl3T > has been neglected. 

Next, including nonlinear terms, we note that the lowest 
order terms (i.e. those proportional to B 2 and v 2 ) do not con- 
tribute directly to the nonlinear coupling, and thus we must 
include nonlinear terms of higher order in both the amplitude 
and the 1 / co c expansion. We thereby obtain 
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The linear relation between the velocity and density perturba- 
tions reads pi = (po/sin0)(9v y /(9^, 
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i.e. Eq. ([TBT l reduces to 
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which is a modified Korteweg de Vries (MKdV) equation with 
a focusing type of nonlinearity. As is wellknown (see e.g. 
Ref. B40I1 ). this equation admits sech-shaped soliton solutions 
where the product of the amplitude and the width is a constant, 
according to 
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Astrophysical environments may exhibit extreme fields. 
Neutron stars have surface magnetic field strengths of the or- 
der of 10 6 — 10 9 T MM], while magnetars field strengths can 
reach 10 10 - 10 11 T [36], coming close to energy densities 
corresponding to the Sch winger limit 112011 : here, in the vicin- 
ity of magnetars, the quantum vacuum becomes fully nonlin- 
ear. However, more moderate, but still very strong, fields ap- 
pear at a distance from the surface of magnetized stars, and 
often in conjunction with a pair plasma due to cascading pro- 
cesses 13511 . Let us therefore consider a specific example of 
solitons in a pulsar environment. We then take the unper- 
turbed magnetic field as Bq ~ 10 T. For a pulsar period of 
P = 1 s, and a multiplicity n = 10 114 111 , the Julian-Goldreich 
expression « JG = 7 x 10 15 (0.1 s/P)(B/10 8 T)m- 3 for the pair 
plasma density gives po ~ 10~ 15 kg/m 3 lf35ll . Furthermore, as 
supported by the work of e.g. Ref. 14211 . we let the pair plasma 
temperature be moderately relativistic, i.e. T ~ QAmc 2 . Fi- 
nally, we let the Alfven waves propagate almost parallell to 
the magnetic field, at an angle 9 ~ 0.2rad, and have a veloc- 
ity amplitude v y o ~ 10 3 m/s. We then find that the width of 
the soliton is d = [(^ B Bo«c) 2 v 2 sin 2 9/3mTc^cos 3 0]" 1 / 2 ~ 
10~ 2 m. We note that variations of the parameters in the ex- 
pression for the soliton width may change the size signif- 
icantly, but this may also require other modifications to be 
made, such as the inclusion of relativistic effects [46]. How- 
ever, at this stage, we conclude that the existence of local- 
ized spin structures in the pair plasmas surrounding pulsars is 
likely, and further generalizations are left for future studies. 



V. SUMMARY AND DISCUSSION 

In the present paper we have studied the dynamics of an 
electron- positron pair plasma in the MHD limit, and included 
the spin properties of the constituents. A closed set of one- 
fluid equations have been obtained, resembling the standard 
MHD equations but including both a standard quantum force 
(from the so called Bohm potential), as well as a number of 
new terms related to the particle spins. The spin terms are of 
particular importance for strongly magnetized plasmas and for 
low temperature plasmas, when the spins are aligned with the 
magnetic field. We stress that the spin terms can have rather 
different properties than the usual terms in the MHD equa- 
tions. As a consequence, it turns out that the spin force can 
be of importance even when its magnitude is smaller than the 
usual j x B- force. In order to demonstrate this property, we 
have studied the special case of weakly nonlinear shear Alfven 
waves in the one-dimensional limit. Due to the spin effects, 
such waves may be governed by a MKDV-equation, leading 
to wave steepening and subsequent soliton formation. By con- 
trast, the nonlinearity cancels to all orders for such waves, if 
the spin terms are omitted. The results of the present paper can 
be of particular significance for astrophysical pair plasmas in 
the vicinity of pulsars and magnetars lU7llT8ll . as well as low- 



temperature laboratory plasmas [43], and nano-structured sys- 
tems Q. 
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